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EQUATIONS 
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Abstract. We study the phase space of the evolution equation 

ht = -{f{h)h xxx ) x - (g(h)h x ) x 

by means of a dissipated energy (a Liapunov function). Here h(x,t) > 0, and at ft, = the 
coefficient functions / > and g can either degenerate to 0, or blow up to oo, or tend to a 
nonzero constant. 

We first show all positive periodic steady states are 'energy unstable' fixed points for the 
evolution (meaning the energy decreases under some zero-mean perturbation) if (g/f)" > 
or if the perturbations are allowed to have period longer than that of the steady state. 

For power law coefficients (f(y) — y n and g(y) = By" 1 for some B > 0) we analytically 
determine the relative energy levels of distinct steady states. For example, with m—n G [1, 2) 
and for suitable choices of the period and mean value, we find three fundamentally different 
steady states. The first is a constant steady state that is nonlinearly stable and is a local 
minimum of the energy. The second is a positive periodic steady state that is linearly 
unstable and has higher energy than the constant steady state; it is a saddle point. The 
third is a periodic collection of 'droplet' (compactly supported) steady states having lower 
energy than either the positive steady state or the constant one. Since the energy must 
decrease along every orbit, these results significantly constrain the dynamics of the evolution 
equation. 

Our results suggest that heteroclinic connections could exist between certain of the steady 
states, for example from the periodic steady state to the droplet one. In a companion article 
we perform numerical simulations to confirm their existence. 



1. Introduction 

We study the evolution equation 
(1) th = -(f(h)h xxx ) x - {g{h)h x ) x . 

This is the one dimensional version of h t = — V ■ (/(/i)VA/i) — V • (g(h)'Vh), which has been 
used to model the dynamics of a thin film of viscous liquid. The air/liquid interface is at 
height z = h(x,y,t) > and the liquid/solid interface is at z = 0. The one dimensional 
equation ([I]) applies if the liquid film is uniform in the y direction. 

The fourth order term in the equation reflects surface tension effects and the second order 
term can reflect gravity, van der Waals interactions, thermocapillary effects or the geometry 
of the solid substrate, for example. Typically f(h) = h 3 + f3h p where < p < 3,(3 > 0, and 
g(h) ~ ±h m as h — > 0, where m G R. In certain applications g(h) changes sign at some 



positive h. We refer to 22, 26 1 for reviews of the physical and modeling literature. 
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The extensively studied Cahn-Hilliard equation ^] also has the form flU), with / = 1 



and g(/i) = 1 — 3h 2 . See [0,0,0, 0] for further references on the Cahn-Hilliard equation. 

Equations like (0) are of mathematical as well as physical interest: for example, Bertozzi 
and Pugh [[| conjectured that blow-up (||/i(-,i)||oo — * oo) is possible in some cases (e.g. if 
f(h) = h n ,g(h) = h m with m > n + 2), and they have proved || finite time blow-up for 
f(h) = h and g(h) = h m when m > 3. In |T7|, §8] we related the steady states and some of 
their properties to this blow-up conjecture. 

Background and goals. In [0 we proved linear stability and instability results for the 
positive periodic steady states of ([[]). Periodicity should not be regarded as a constraint, 
since if /, g > then positive bounded steady states must be periodic or constant, by fPT , 



Theorem B.l]. And periodic steady states do exist for many equations of type ([I]), by the 
methods of |T7|, §2.2] or JIB], ^5], for example. Compactly supported 'droplet' steady states 
only exist, though, if g/ f satisfies additional constraints |17|, §2.2], and can have relatively 
low regularity at the contact line. 

In this paper we concentrate mostly on positive periodic steady states, and on droplet 
steady states with zero contact angles. Our main investigative tool is the energy 



s(h(.,t))= r 

Jo 



h x (x,t) 2 - H(h(x,t)) 



dx: 



here h(x,t) is a smooth solution of (P that is X-periodic in x, and H(y) satisfies H" = g/f. 
This energy is strictly dissipated: (d/dt)£(h(-,t)) < with equality if and only if h is a 
steady state (see § |2~Tf) . Thus the energy is a Liapunov function for the evolution. 
We address two questions about the energy landscape of the evolution (p. 

1. Which steady states are local minima of the energy? Which are saddle points? 

2. Among steady states having the same period and the same area (i.e. fluid volume), 
which has the lowest energy? 

Answering these questions will help clarify the phase portrait of the evolution. For example, 
if all small zero-mean perturbations of a steady state can be shown to raise the energy, then 
the steady state might be asymptotically stable: it might be that every smooth solution 
starting from near the steady state relaxes back to the steady state as t — > oo. But if some 
zero-mean perturbation decreases the energy, then asymptotic stability definitely fails. 

Our requirement that the perturbations have zero mean seems reasonable from a physical 
standpoint, because it corresponds to a disturbance of the fluid that alters the profile without 
adding additional fluid. Mathematically it is reasonable because the evolution equation (0) 
preserves volume for spatially periodic solutions: J h(x,t)dx = J h(x,0)dx for all time t. 
Thus zero-mean perturbations allow the possibility of relaxation back to the original steady 
state, while nonzero-mean perturbations do not. 





Figure 1. Four types of steady state. 
A sketch of definitions and results. Take X > 0. If there is an X-periodic zero-mean 
perturbation v such that £(h SB + ev) < £(h ss ) for all small e > 0, then we call the steady 
state h ss 'energy unstable' at period X. If instead £(h ss + ev) > £(h ss ) for all sufficiently 
small e > 0, for each X-periodic zero-mean perturbation v, then we call the steady state 
'energy stable' at period X. (Some authors call this formal stability |16|].) It is conceivable 
that a steady state might be energy stable and yet not be a local minimum of the energy. 
Our main stability results, in Section |2|, are roughly stated as follows. 

• Theorem |]. For positive periodic steady states, linear instability implies energy insta- 



bility. Hence our linear instability results in |18| imply that every positive X*-periodic 
steady state is energy unstable at periods X = 2X*, 3X*, . . . , and is also energy unsta- 
ble at period X* if g/f is a strictly convex function. 

• Theorems Q-|3|. Further, for the 'power law' coefficients f(y) = y n and g(y) = By m with 
B > 0, we completely characterize energy stability at period X* even when g/ f is not 
convex, that is, when m — n G (0, 1). 

In Section [2.6| we explain how these results for periodic boundary conditions cover the case 
of Neumann ('no flux') boundary conditions as well. 

Then in Section |] we determine the relative energy levels of three different kinds of steady 
state: constant steady states, positive periodic steady states, and zero contact angle droplet 
steady states. Figure || illustrates these three steady states, as well as showing a nonzero 
contact angle droplet steady state (about which we say little in this paper). 

We have found it too difficult to determine energy levels when working with arbitrary 
coefficients / and g, but have obtained fairly complete answers for power law coefficients. 
This provides at least some insight into the general case. Further insight comes from the 
work of Grinfeld and Novick-Cohen [14| on the energy levels of steady states for the Cahn- 
Hilliard equation, which has non-power-law coefficients. The earlier work of Mischaikow 
applies to a variety of gradient-like bistable equations. 
Our main energy level results, in Section |^, for f(y) = y n and g(y) = By m are: 

• Theorem || If m < n or m > n+1 then a positive periodic steady state always has higher 
energy than the constant steady state with the same mean value. For n < m < n + 0.75 
our analytical and numerical work suggest the positive periodic steady state has lower 
energy than the constant steady state. 



• Theorem When m ~ n + 0.77 there can be two steady states with the same period 
and area, h BB i and h SS 2, with min x h BB \(x) < min x h SS 2(x); we essentially prove h SS 2 is 
energy unstable and has higher energy than h ss i, which is energy stable. 

• Theorem 0. If m — n G (—2, 0) U [1, 2) then a positive periodic steady state always has 
higher energy than a zero contact angle droplet steady state with the same mean value. 

• Theorem |TI]. The constant steady state can have higher energy than the zero contact 
angle droplet steady state with the same mean value. When m — n G [1,2), for example, 
a mountain pass scenario can occur, in which the constant steady state is a local 
minimum of the energy, the positive periodic steady state is an energy unstable saddle, 
and there is a zero contact angle droplet having lower energy than either of them. 

For example, Theorem ^ with m = —1 < n = 3 covers the 'van der Waals' case 
(2) h t = -(h 3 h xxx ) x - A{h- 1 h x ) x 



with A > 0. This equation has been studied by a number of other authors, e.g. j?J, 30], |3T, [33 
mostly with regard to similarity solutions and film rupture (where the solution goes to zero 
in finite time, at some point). Our numerical work on this equation pays particular 
attention to behavior near steady states, and how this evolves into rupture. 

Our energy level results suggest possible basins of attraction around the stable steady 
states, and possible heteroclinic connections between steady states. In the companion arti- 
cle 0, §4], we investigate such possibilities with numerical simulations. For example, when 
m ~ n + 0.77 as in Theorem ^, we find robust heteroclinic connections between the unstable 
positive periodic steady state h SS 2 and the stable one h BS i. For the mountain pass scenario 
in Theorem [U], we find that perturbing the saddle point (the periodic steady state) in one 
direction leads to relaxation to the constant steady state and perturbing in the opposite 
direction gives apparent relaxation to a droplet. A similar dichotomy was found for ax- 
isymmetric surface diffusion by Bernoff, Bertozzi and Witelski, [||, p. 744], with perturbed 
unduloids either relaxing to a cylinder or else pinching off in finite time. 

We also present in [IS, §5] simulations suggesting that small changes in the 'mobility' 



coefficient / do not break heteroclinic orbits, but can affect whether or not the solution 
remains positive as it evolves. 

We discuss some of these conclusions and future directions further in Section ||. 

Terminology. We write Tx for a circle of circumference X > 0. As usual, one identifies 
functions on Tx with functions on R that are A-periodic and calls them even or odd according 
to whether they are even or odd on R. 

A positive periodic steady state is assumed to satisfy the steady state equation classically. 
A droplet steady state h ss (x) (see Figures |l|c and d) is by definition positive on some interval 
(a, b) and zero elsewhere, with h ss G C^a, &]; we require h ss to satisfy the steady state 
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equation on the open interval (a, b) only, and to have equal acute contact angles: < 
h' ss (a) = —h' ss (b) < oo. (Throughout the paper, if a function has only one independent 
variable then we use ' to denote differentiation with respect to that variable: h' ss = (h S8 ) x .) 

We say a droplet steady state h ss has 'zero contact angle' if = h' ss (a) = —h' ss (b), and 
'nonzero contact angle' otherwise. A 'configuration' of droplet steady states is a collection 
of steady droplets whose supports are disjoint. For more on the steady states and their 
properties, see [17[ . 



2. Energy stability for periodic steady states 

We assume throughout this section that f(y) and g(y) are C 2 -smooth for y > 0, and that 
/ > 0. Define 

9 

r = T 

Take A > 0. 

We investigate stability and dynamical questions by means of a Liapunov energy. A few 
of the theorems follow directly from our linear stability results in [^8j, but most are quite 
different and complementary. 

2.1. Definition of the energy, and of energy instability. The energy function for the 
evolution equation ([I]) is defined for I e H l (Tx) to be 



(3) £(£) 



o 



dx, 



where H(y) is a function with H" = r = gj f. 

To verify the energy £ is a Liapunov function for the evolution (|l|), suppose h(x,t) is a 
positive smooth solution of (0) that is A-periodic in x. Bertozzi and Pugh || §2] observed 



(generalizing JIT], [12], |27[) that £ is dissipated by the evolution: 

j t £(h(-,t)) = - ^ [f(h(x, t))h xxx (x, t) + g(h(x, t))h x (x, t)f dx < 0. 

The dissipation is strict at each time t unless f(h)h xxx + g(h)h x = for all x. For smooth 
positive periodic solutions this occurs only when h(-,t) is a steady state. 

Let h ss G C 4 (T X ) be a positive periodic steady state of (|l|). It is easy to see (cf. formula 
(fL5D) that h ss is a critical point for the energy £, with respect to zero- mean perturbations. 

Definition. Call h ss an energy unstable critical point (with respect to zero-mean perturba- 
tions at period A) if there exists a smooth A-periodic perturbation u(x) with mean value 
zero such that 

£(h ss + eu) < £{h ss ) for all small e > 0. 



That is, small perturbations in the direction u decrease the energy. (Some authors call this 
forma/instability Q.) 

An energy unstable steady state is necessarily a saddle point in the energy landscape, 
since £ is increased by the perturbation u(x) = e cos(2iikx/X) for k ;§> 1. 

Energy unstable steady states are not asymptotically stable in H 1 (Tx), in the following 
sense: suppose h(x,t) is a positive smooth solution with initial data h ss + eu; then h(-,t) -f+ 
h ss (-) in if^Tx) since for all t, £(h(-,t)) < £(h(-,0)) = £(h ss + eu) < £(h ss ) (convergence 
in H 1 would imply convergence in L°° and hence convergence of the energy). In fact, h(-,t) 
cannot converge to any translate of h ss , for the same reason. 

The last paragraph extends to nonnegative weak solutions if they also dissipate the energy. 
Weak solutions must sometimes be considered because solutions of (P that are initially 
positive might not always remain so, and where they go to zero they can lose regularity. 
See Bertozzi and Pugh || |6| for existence of nonnegative weak solutions that dissipate the 
energy. 

2.2. Energy instability results. In [fT8| , §2] we linearized the evolution equation ([I]) around 
the positive periodic steady state h ss and then reduced the linear stability question to deter- 
mining the sign of the first eigenvalue of a certain self-adjoint fourth order linear operator. 
We will not repeat the linearization here, or re-state the linear stability results of [18|. How- 
ever we warn readers that when we say a steady state is 'linearly stable', we are including 
the neutrally stable case in which the first eigenvalue of the linearized operator is zero. This 
is unavoidable: the operator always has a zero eigenvalue in its spectrum, corresponding to 
an infinitesimal translation of the steady state in space (the evolution equation is translation 
invariant). 

The next theorem states that if a steady state is linearly unstable then it is energy unstable. 
Also, we present some unstable directions, when g/f is strongly convex. 

Theorem 1. Let f,g G C 2 (0, oo) with f > 0. Take X > and suppose h ss e C 4 (IR) is an 
X -periodic positive steady state of (]![). 

If h ss is linearly unstable with respect to zero-mean perturbations at period X, then it is 
also energy unstable at period X . 

In particular, h ss is energy unstable at period X if it is non-constant and either: the least 
period of h ss is X/j for some integer j > 2 or else r = g/f is convex (r" > 0) and non- 
constant on the range of h ss . For example, if h ss is non-constant and r" > 0, then h ss is 
energy unstable in the directions u = ±h' ss and i/i^. 

This is proved in Section |4.1| . Note for example that the theorem covers the van der Waals 
evolution (|J), since there r(y) = Ay~ l /y z = Ay~~ 4 is strongly convex. Thus all positive 
periodic steady states of the van der Waals evolution are energy unstable. 



These energy unstable steady states (which we observed above are not asymptotically 
stable) are presumably nonlinearly unstable, in general, but we cannot prove this. 

In Theorem H] we have assumed f,gE C 2 (0, oo), which is generally the case for the thin film 
equations that are our main motivation. But our arguments are all local (involving only small 
perturbations of the steady state), and so the theorem still holds if the coefficient functions 
f(y) and g(y) are defined and C 2 merely for ^-values in a neighborhood of [h ssmin , h ssnirLrl 



2.3. Review of power law steady states and their rescalings. We now turn to power 
law coefficients: f(y) = y n and g(y) = By m for some exponents n, m G R and some positive 
constant B > 0. Here 

r(y) = By*' 1 

where 



m — n + 1 



This exponent q determines many properties of the steady states, including (usually) their 
linear stability. 

The evolution equation ([[]) becomes 

ht = —{h n h xxx ) x — B(h m h x ) x . 

To state our results on energy stability for this power law evolution, we first review some 
properties of the steady states and explain how to rescale them to solutions of a canonical 
nonlinear oscillator ODE, given as equation ([8|-^]) below. 

We start with a non-constant positive periodic steady state h ss e C A (Tx) of the general 
evolution ([!]). The steady state condition for (0) integrates to give f(h BS )h' s " + g(h ss )h' ss = C 
for some constant C . The least period of h ss is X/j for some integer j ' > 1. 

One finds that the constant C (the flux) equals zero, by dividing f(h ss )h'^ + g(h ss )h' ss = C 
by f(h ss ) > and integrating over a period (cf. |^7|, p8|). Hence the steady state satisfies 



(4) hZ + r(h ss )h' ss = 0. 



[If h ss were a droplet steady state then again C = 0, by |T7|, Theorem 2.1], and equation (|) 
would hold wherever h ss is positive.] 

Integrating, the steady states have a nonlinear oscillator formulation: 

(5) hi + H'(h BS ) = 

holds wherever h ss is positive. Here H(y) is a function with H" = r = g/f; if we regard x 
as a 'time' variable then ^h' ss (x) 2 + H(h ss (x)) is a conserved quantity. 



Returning to the power law evolution, remember r(y) = By q 1 . Thus for q ^ we can 
write the steady state equation @ as 

Bh q - D 
Q 

for some constant D. For q = the analogous equation is /i" s + Blogh ss — D = 0. This 
oscillator equation involves three constants: q, B, and .D. We remove B and -D by rescaling: 
let 



(7) k(x) 



(f) 1/9 ^((f) 1/29 ^) ; **o, 

e-^ ss (e B / 2 ^), g = 0. 

For g 7^ this rescaling uses that -D > 0, by 0, §3.1]. [A different rescaling would be used 
to study droplet steady states fl7|, §§3.2,4].] The steady state equation (|J) rescales to 

(8) k" + ^—^ = 0, q^O, 

(9) k" + log k = 0, q = 0. 

Differentiating, we find for all q that k'" + k q ~ l k' = 0, and so k satisfies (k n k"' + k m k')' = 0, 
i.e. it is a steady state of k t = —ik n k xxx ) x — (k m k x ) x . 

Since h ss is non-constant, positive and periodic, we know k"(xo) > for some point Xq. 
Evaluating (^-0) at xq shows the minimum value of k is less than 1. Also A;'(0) = since 
(after a suitable translation) h ss has its minimum at x = 0. Introducing the notation k a for 
the solution k that has minimum value a G (0, 1), at x — 0, we have 

(10) < k a (0) =a<l, £4(0) = 0. 

Thus every steady state h ss can be rescaled to a fc a , as above. Conversely, for each q 6 R 
and a G (0, 1) there exists a unique smooth positive periodic k a satisfying equations (||-f5D 
and ( |I0] ) (see |17|, Proposition 3.1]). The same holds for a = when g > —1, except that 
ko may be only C 1 -smooth at x = (see |T7|, Theorem 3.2]). To illustrate, Figure ^ plots 
the steady states k a over two periods, for q = 3 and eight a- values between and 1; see ]T7|, 
§6.1] for details. 

Note that the map (a, x) \— > k a (x) is C°°-smooth for (a, x) £ (0,1) x IR, by an ODE 
theorem giving smooth dependence on the initial data [[TJ|, Ch. V §4]. We write 

P = P(a) and A = A(a) 

for the least period of k a and for the area under its graph, A = ' k a (x) dx, respectively. 
Then P and A are smooth functions of a that approach 2tt as a — > 1, by Lemma 6]. 
The function 



E(a) := Piaf^Aiay- 1 = P(a) 2 [A(a)/P(a)] 



q-l 
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FIGURE 2. Steady states k a (x), when q = 3. 

determines whether the steady state is energy unstable or stable, in several results below. 

The above rescaling ideas are a useful tool throughout the paper. We hope this tool does 
not obscure the fact that stability and energy level properties for equations of type (|1|) seem 
to be determined by the period map of a family of steady states h SB with fixed area but 
varying amplitudes and periods, or alternatively the area map of a family of steady states 
with fixed period but varying amplitudes and areas. This is how one should think of the 
function E(a); see [|18| , §6.3] for more on this. The same underlying idea appears in the work 
of Grinfeld and Novick-Cohen [[L4| on the Cahn-Hilliard equation, an equation which is not 
amenable to rescaling in the same way. 

2.4. Energy in/stability for the power law evolution. 

Theorem 2. Let h ss G C 4 (M) be a non- constant positive X -periodic steady state of the power 
law evolution ht = —{h n h xxx ) x — B(h m h x ) x . Translate h ss to put its minimum at x = so 



that h ss rescales to k a for some a G (0, 1), as in §2.5. 

If q < 1 or q > 2 then h ss is energy unstable in the directions u = ih' 8s and u = ±h" s 
If q = 2, or if q > 1 and E'(a) > 0, then h ss is energy unstable. 



We prove the theorem in Section |4.2| . Its first statement follows immediately from The- 
orem [I], since r(y) = By q ~ l is strongly convex (r" > 0) when q < 1 or q > 2. The final 
statement of the theorem certainly applies when q > 2, since then E' > by |]18, Theo- 



rem 11]. Computational studies (J7[, §6.1] suggest E'(a) > for all a G (0, 1) except when 
q G [1, 1.794] (approximately); Figure || plots E(a) for certain of these g-values. 

To help explain the appearance of the criterion E'(a) > for energy instability, in the 
theorem, we refer the reader to 0, §6.3] and to j0| §3] for a 'bifurcation diagram' interpre- 
tation of the function a i— > E(a) in terms of a family of steady states h BS with fixed area but 
varying amplitudes and periods, or fixed period but varying amplitudes and areas. 

The case 1 < q < 2 with E'(a) < is not covered by Theorem ^. By Jl8|, Theorem 9] we 
know h ss is linearly stable in this case. We have been unable to prove nonlinear stability, but 
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we do prove in Theorem |3| below that if E'(a) < then small perturbations of h ss in every 
possible direction strictly increase the energy; this is consistent with nonlinear stability. 

Definition. Let h SB G C 4 (Tx) be a non-constant positive periodic steady state of ([!]). 
Call h SB energy stable (with respect to zero-mean perturbations at period X) if for each 
u G H 1 (Tx) \ {0} with u dx = we have 

£(h ss + eu) > £{h ss ) for all small e > 0. 

(Some authors would call this formal stability p~6|] .) 

A steady state might perhaps be energy stable without being a local minimum of the 
energy. (For example, the function f(x,y) = (y — x 2 )(y — 3x 2 ) on M 2 has the property 
that origin is a local minimum point on each straight line through the origin, though it is 
not a local minimum in R 2 .) Another cautionary note is that the energy is insensitive to 
translation: in particular, a steady state and its translates have the same energy, consistent 
with the translation invariance of the evolution equation itself. Numerical simulations in our 



article [19, §4.4] demonstrate that perturbations of an energy stable steady state can evolve 
towards a translate of that steady state. This suggests that any asymptotic or nonlinear 
stability result that can be proved will hold only up to translation. 



In Section 4.3 we prove: 



Theorem 3. Let h ss G C 4 (1R) be a non- constant positive periodic steady state of the power 
law equation h t = — (h n h xxx ) x — B(h m h x ) x , with h ss having least period X . Translate h ss to 
put its minimum at x = so that h ss rescales to k a for some a G (0, 1), as in §[0. 
If 1 < q < 2 and E'(a) < then h ss is energy stable. 

The hypothesis E'(a) < seems numerically to hold for all a G (0, 1) when 1 < q < 1.75, 
as indicated by Figure 3 in |ffS| . 



We have no energy stability result when 1 < q < 2 and E'(a) = 0; fortunately it seems 



E'(a) = for at most one a-value, for each q, as shown numerically by Figures 3-5 in [|T 
Theorem 9 of does imply linear stability when E'(a) = 0, but [18|, Theorem 10(b)] shows 
the space of neutrally stable directions is two dimensional (rather than one dimensional as 
when E\a) < 0) and this might perhaps lead to instability. 
We next address the q = 1 case. 

Lemma 4. Let q = 1 (i.e. m = n) and suppose h ss G C 4 (R) is a non-constant positive 
periodic steady state of h t = —{h n h xxx ) x — B(h n h x ) x with least period X, and translate so 
that h ss has its minimum at x = 0. Then h ss is not asymptotically stable in if 1 (Tx) with 
respect to even perturbations. 
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Proof. The steady state equation (||) with q = 1 has general solution h ss (x) = D/B + 
ccos(v^Ba;), where we have used that h ss has an extremum at x — 0. Hence the period is 
X = 2tc/^/B, and for all small e we see that the perturbed function h ss (x) + e cos(VBx) 
is another positive periodic steady state solution. Thus h SB is not asymptotically stable in 
H 1 (T X ) with respect to even perturbations. □ 

The steady state is of course not asymptotically stable with respect to general perturba- 
tions, since one can always perturb by translating the steady state a small distance (this 
remark applies for all q). However we are really not interested in such simple transla- 
tional motion. Also, translational perturbations are not even permitted under the Neumann 
boundary conditions that we consider below, in Section |2.6| . 

To summarize, when q = 1 the positive periodic steady states are linearly neutrally stable 
with respect to zero-mean perturbations of the same period, by |TJ[ or [jTSL Lemma 8], and 
are not asymptotically stable with respect to 'even' perturbations, by the above lemma. 
Our numerical simulations in [19|, §4.3] and those of ]12| suggest that a wide range of small 



perturbations yield solutions relaxing to nearby positive periodic steady states, suggesting 
they are nonlinearly stable. 

In the companion paper §4] we illustrate the above stability and instability theorems 
for the power law evolution with a variety of numerical simulations. There we find not only 
the short time behavior suggested by the energy (in) stability results, but also some longer 
time limits that are suggested by the energy level results in Section |3|. 

2.5. Odd perturbations. Returning momentarily to general coefficients / and g, in Section 



4.4 we prove the energy increases under odd perturbations, when r is concave. 



Theorem 5. Let h ss £ C A (R) be a non-constant positive periodic steady state of with 
coefficient functions f, g £ C 2 (0, oo),/ > 0. Suppose h ss has least period X, and translate 
h ss to put its minimum at x = 0. 

If r = g/f is strongly concave (r" < 0) then for every nontrivial u £ H 1 (Tx) that is odd 
and is such that h ss + u > 0, we have S(h ss + u) > £(h ss ). 

The theorem is global since the perturbations are not required to be small, and it is 
consistent with asymptotic stability (although convergence to a translate of h ss seems more 
likely than convergence to h ss itself). 

Theorem [| applies in the power law case with 1 < q < 2,r(y) = y q ~ l . Another example 
with strongly concave r is an equation eqn. (24)] describing the dynamics of a population 
of aphids, for which f(y) = y, g(y) = y — c and r(y) = 1 — c/y. 

2.6. Relation between the periodic and Neumann stability problems. Suppose h ss 
is an even A-periodic steady state of the evolution equation ([!]) with extrema at x — 
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0, ±X/2, so that h' ss = h"l = at these points. As we observed at the end of |18|) 
§2.5], linear instability of h ss with respect to periodic boundary conditions on (—X/2,X/2) 
is equivalent to linear instability with respect to Neumann ('no flux') conditions on the 
half-interval (0,X/2); these Neumann conditions are: h x = h xxx = at x = 0,X/2. 

The energy of a positive smooth solution is still dissipated in the case of Neumann bound- 
ary conditions, and obviously energy instability of h ss in an even direction u(x) for the 
periodic problem on (— X/2, X/2) is equivalent to energy instability in the direction u(x) for 
the Neumann problem on (0,X/2). (If the perturbation u(x) is even and has mean value 
zero on (—X/2, X/2) then it has mean value zero on (0,X/2) as well.) Similarly, 'periodic' 
energy stability in all even directions on (—X/2, X/2) is equivalent to 'Neumann' energy 
stability in all directions on (0,X/2). 

In short, for the Neumann problem on (0,X/2), the stability result in Theorem |3| still 
holds, and the instability claims involving ±/i" s in Theorems [l] and also still hold, since 
these are even functions; the claims involving ±h' ss do not carry over, since those are odd. 



3. Relative energy levels of periodic, constant and droplet steady states 

In this section we investigate the phase space of the power law equation h t = —(h n h xxx ) x — 
B(h m h x ) x by comparing the value of the energy at the positive periodic, constant and zero- 
angle droplet steady states. Let X > and recall q — m — n + 1. 

3.1. Positive periodic vs. constant steady states. The fluid volume h(x,t) dx is 
conserved by the evolution, under periodic boundary conditions, and so the mean value 

i r x 

h := — / hdx 



XJ 

is constant in time. Suppose the initial data h(-, 0) arises from a small zero-mean perturbation 
of h ss . It is natural to ask whether h can converge (while staying positive and smooth) 
towards the constant steady state h ss , as t — > oo. This cannot happen if S(h ss ) < S(h ss ) and 
also h(-, 0) is close enough to h ss so that S(h(-, 0)) < £(h ss ), because the energy is dissipated 
by the evolution. 

Theorem 6. Let h ss G C 4 (1R) be a non-constant positive periodic steady state of the power 
law equation h t = —(h n h xxx ) x — B(h m h x ) x , with least period X. Translate h ss to put its 



minimum at x = so that h ss rescales to k a for some a € (0, 1), as in §2.5 
If q>2 or q < 1 then £(h ss ) > £(h ss ). 
Ifq = l then £(h ss ) = £(h ss ). 

Ifl<q<2 and E' > on (a, 1) then £(h ss ) > E(h ss ). 
Ifl<q<2 and E' < on (a, 1) then £(h ss ) < E(h ss ). 



-3 it was observed numerically in 
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}3] 



The theorem is proved in Section |5.1| . For q 
that S(h BB ) > S(h ss ). 

When 1 < q < 1.75, numerical evidence in Figure 3 of |1| suggests E'(a) < for all 
a. If this is true, then S(h ss ) < S(h ss ) by the theorem and so there can be no heteroclinic 



connection from h ss to the constant steady state h ss ; indeed in |T9, §4.4] we find numerically 
for q = 1.5 that h ss is asymptotically stable, up to translation. 

For 1.795 < q < 2, Figure 5 of @ suggests E'(a) > for all a and so £(h ss ) > £(h ss ) 
by the theorem. Thus when q > 1.795 or q < 1, the instability result Theorem ^ and the 
energy level result Theorem [I] together lead us to suspect the existence of a heteroclinic 
connection from h ss to h ss . In Figured we present numerical simulations of such heteroclinic 





Figure 3. q = 2.5 and n 
dashed line: initial data. 



1. 



Figure 4. q = 1.768, n 
dashed line: initial data. 



orbits, taken from the companion article [T^, §4.6]. The top part of the figure presents an 
orbit connecting the positive periodic steady state to the constant steady state. The bottom 
part relates to § |57^ and presents a solution connecting a perturbation of the same positive 
periodic steady state in finite time to a droplet profile (which will not in general be a steady 
state). 

In particular the van der Waals equation (§) has q = — 1 — 3 + 1 = —3, and so by the first 
claim of Theorem |6], the energy of any positive periodic steady state is greater than that of 
the constant steady state. This was observed numerically by Witelski and Bernoff |32], §3]. 
See |32| and our companion paper [19|, §4.1] for numerical simulations of this equation 
and a discussion of the literature. 

For q G (1.75, 1.794], see the remarks around Theorem ^. 

For the Cahn-Hilliard equation, the analogue of Theorem |6] (comparing energy levels of 
non-constant and constant steady states) can be found in Grinfeld and Novick-Cohen's work 
14] , Theorem 4.1]; further, |1^, §7] discusses a number of results on existence of heteroclinic 

3.4]. 



connections. See also |21 



3.2. Positive periodic vs. droplet steady states. We do not yet have stability results 
for droplet steady states. But here we show under certain conditions that the energy of a 
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zero-angle droplet steady state must be lower than that of a positive periodic steady state h ss 
whose period exceeds the length of the droplet. In our article |19|, §4] we show numerically 
that in these cases, the droplet seems to be strongly attracting, with heteroclinic orbits from 
h ss towards the droplet steady state as shown in the bottom part of Figure [| 

Theorem 7. Let h ss G C 4 (M) be a non- constant positive periodic steady state of the power 
law equation h t = —(h n h xxx ) x — B(h m h x ) x , with h ss having least period X and area A ss = 
Jq h ss dx. 

If -1 < q < 1 or 2 < q < 3, or if 1 < q < 2 and E'(a) > for all a G (0, 1), then there 
exists a zero contact angle droplet steady state h ss with length X < X and area A = A ss . 
Furthermore S(h ss ) < S(h ss ). 



We prove the theorem in Section p?2 
Remarks 



1. Steady states with zero contact angle can occur only for q > —1, by [T7|, §2.2]. Hence 
we do not consider q < — 1 in the theorem; indeed variational methods show there can be 
droplet steady states that do not have acute contact angles p4| . 

2. The theorem presumably applies to 1.795 < q < 2, since it seems numerically that 
E' > for those g-values. For q > 3 we think a similar theorem might hold but with h ss 
being a configuration of disjoint zero-angle droplets. 

3. Since £(h ss ) > £(h ss ), there might be an orbit from h ss to h ss . This orbit might 
describe a positive solution that converges to the nonnegative droplet profile as t — > oo, or 
it might describe a positive solution that loses positivity in finite time and then approaches 
the droplet as a nonnegative weak solution. For — 1 < q < 1 and 2 < q < 3, Theorems |6] 
and are consistent with the idea that the unstable positive periodic steady state h ss and 
its stable manifold form a separatrix between the basin of attraction of the constant steady 
state and the basin of attraction of the droplet steady state. In particular, after perturbing 
h ss in one direction one seems to find a solution that converges to the constant steady state, 
while perturbing in the opposite direction often yields a solution that converges to a droplet 
profile. We present some numerical evidence for such behavior in Figure and discuss 
this at length in the companion article |19|, §4]. See also the 'mountain pass' remark after 
Theorem [11]. 

The preceding theorem and its Remarks are definitely not valid for q G (1,1.75], since 
for these q- values we show there does not even exist a zero contact angle steady state with 
length less than X: 

Theorem 8. Let 1 < q < 1.75 and suppose E'(a) < for all a G (0,1). Let h ss G 
C 4 (Tx) be a non-constant positive periodic steady state of the power law equation h t = 
— (h n h xxx ) x — B(h m h x ) x , with h ss having least period X and area A ss = h ss dx. Assume 
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h ss is nonnegative and piecewise-C 1 on Tx, has area A ss , and is smooth on the set where it 
is positive and satisfies there the 'nonlinear oscillator' steady state equation (Qj. 

Then h ss is either constant or is a translate ofh ss , or is a configuration of nonzero contact 
angle droplet steady states. Specifically, h ss cannot be a zero contact angle droplet steady 
state. 



We prove the theorem in Section |5]3|. Note that the hypothesis E' < seems to hold for 
1 < q < 1.75, by the numerical evidence in Figure 3 of fl8|j . 

Finally, for q G (1-75, 1.794] approximately, we know by the analytical and numerical 
work in |I7| , §5.1] that there can be two positive periodic steady states with the same period 
and area. The next theorem shows that the steady state with smaller minimum value (and 
larger amplitude) is energy stable, and has lower energy than the other one, which is energy 
unstable. This leads us to suspect there exists an orbit from the energy unstable steady 
state to the energy stable one, at least when the steady states have been chosen to have their 
minima at the same location. The bottom plot of Figure § presents a numerical simulation of 
such an orbit, taken from fl9| , §4.5]. The top plot of the figure presents an orbit connecting 
the unstable positive periodic steady state to the constant steady state, similar to the top 
plot of Figure |3|. 

Theorem 9. Assume 1 < q < 2 and there exists a cr u G (0, 1) with E'(a) < on (0, a cr u) 
andE'(a) > on (a cr i t , 1), and assume a i— > aP(a) 2 ^ q ~^ is strictly increasing for a G (0, 1). 

Suppose h ss \ and h SS 2 are non-constant positive periodic steady states of the power law 
equation h t = —(h n h xxx ) x — B(h m h x ) x , with h ss i and h SS 2 having the same least period X and 
same area h ss i dx = h ss2 dx. 

If h ss i(x) has lower minimum value than h SS 2(x), then h ss i is energy stable, h ss2 is energy 
unstable, and S(h ss i) < £(h SS 2). Furthermore, £(h SS 2) > £(h SS 2). 



We prove this in Section |5.4j . The hypothesis about E(a) being first strictly decreasing 
and then strictly increasing is confirmed numerically for q in the interval (1.750, 1.794] by 
1% §6.1]; see Figure |5|. Numerical work also confirms the hypothesis about «P(a) 2 ^ ?_1 ' 



being a strictly increasing function of a for all q G (1, 2). 

Theorem P is analogous to |14], Theorem 4.1(v)] for the transitional and metastable cases 



of the Cahn-Hilliard equation, where Grinfeld and Novick-Cohen show the energy of a 
monotonic 'interface' solution is less than that of a monotonic 'spike' solution having the 
same length and area. 

3.3. Constant steady states: stability results. For any number h > the constant 
function h ss = h is a steady state of the general evolution equation (JJ). We now develop 
analogues for this constant steady state of our earlier stability results. Then in the following 
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0.25 0.5 0.75 1 



Figure 5. Plots of E(a) for q = 1.75,1.76,1.768,1.78,1.79. Top curve: 
q = 1.75, the curves move down as q increases. 



subsection we can compare the energies of the constant and droplet steady states, for power 
law coefficients. 



We start by recalling a characterization of linear instability for steady states. From |18 
Appendix A], an X-periodic positive steady state h BB G C 4 (1R) of (JT]) is linearly unstable with 
respect to zero-mean perturbations at period X if and only if Ti(h BB ) < 0, where 



f X l(u') 2 -r(h ss )u 2 }dx 



11) n(h BB ) := min 



J* u 2 dx 

the minimum here is taken over nonzero u G H 1 (Tx) with udx = 0. The trial function 
u corresponds to a perturbation of h ss . 

Theorem 10. Let f,g G C 2 (0, oo) with f > 0, and write r — g/ f. Let h, X > 0. Then for 
the constant steady state h ss = h of (^j, the eigenvalue T\ in is Ti(h) = (2tt/X) 2 — r(h) . 
The Ti(h)-eigenspace is spanned by sin(27nc/X) and cos(27ra;/X). 

Thus with respect to zero-mean perturbations at period X , the constant steady state is 

linearly unstable if r(h)X 2 > An 2 , 
linearly neutrally stable if r(h)X 2 = 4-7T 2 , 
linearly asymptotically stable if r(h)X 2 < An 2 . 

(a) If r(h)X 2 > An 2 , or ifr(h)X 2 = An 2 and r"(h) > 0, then the constant steady state is 
energy unstable in the directions ± sin(27ra:/A) and ± cos(27rx/A) . 

(b) If r(h)X 2 < An 2 , or if r(h)X 2 = An 2 and r"(h) < 0, then the constant steady state is 
energy stable with respect to zero-mean perturbations of period X . In fact, ifr(h)X 2 < An 2 
then the constant steady state is a strict local minimum of the energy with respect to zero- 
mean perturbations in if 1 (Tx), and h is nonlinearly stable under the evolution ^), in an 
H l (Tx) -sense made precise in the proof. 



We prove the theorem in Section |5.5| . Its linear stability assertions are well-known and are 
included for the sake of completeness. 
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Goldstein, Pesci and Shelley fl2| , §IIIB] used the energy to prove nonlinear instability of 
the constant steady state for the q — 1 case (f(y) = y n ,g(y) = By n ,X = 2n) with either 
2 < B < 4 or B = j 2 for some integer j > 2. 

3.4. Constant vs. droplet steady states. Consider power law coefficients, so that r(y) = 
By q ~ 1 ; then the previous theorem shows that the stability of the constant steady state h is 
determined by whether the quantity Bh Q 1 X 2 is >, =, < 47r 2 . 

Fix X > 0. Does a zero-angle droplet steady state exist with length at most X and with 
the same area hX as the constant steady state? If such a droplet steady state exists, can it 
have lower energy than the constant steady state? 

In this direction, in Section |5.6| we prove: 



Theorem 11. Let h, X > 0, and consider the constant steady state h ss = h of the power 
law evolution equation h t = — {h n h xxx ) x — B{h m h x ) x . 

(a) Suppose — 1 < q < 3. Then there exists a zero contact angle droplet steady state h ss of 
length X < X and area hX if and only if 



(12) 



Bh Q l X 2 > £7(0) =: E {q). 



If such a droplet steady state exists, then £{h) > £(h ss ) if and only if 

(3-«)/« 

= :L(q) (for-l<q<3,q^0) 



(13) Bh q V > A(0) 2 



3 + q 



L(3- 9 )(g + l) 



^ — i 



or Bh X 2 > 4e 2 7r/3 =: L(0) (for q = 0). 

(b) For q = 3, such a steady state h ss exists if and only if Bh X 2 = -£7(0). For q > 3, h 
exists if and only if Bh 9 1 X 2 < -£7(0). For all q > 3, if h ss exists then £(h) < £(h ss ). 

To understand conditions ( |12"D and fllcf) see the plots of E (q) and L(q) in Figures || and |7| 
(constructed using the formulas for A(0), -P(O), -£7(0) in [|T|], §3.1.2]). The graphs of _£7 and L 






Figure 6. Plots of -£7 (g) and L{q). 



Figure 7. Close-up view. 



intersect at q = —1, 1, 3 and at q** ~ 1.775. For 1 < q < q** the figure suggests _£7 > L, and 
so if the droplet steady state in Theorem |TT|(a) exists then it certainly has lower energy than 
the constant steady state. On the other hand, it appears that L > -£7 when — 1 < q < 1 and 
when q** < q < 3; in these cases the energy condition ( |13"D may or may not hold when the 
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existence condition fll2|) holds, so that the energy of the droplet steady state (if it exists) 
might be higher or lower than that of the constant steady state. 

The dashed line at height 4tc 2 in Figures |6] and |7| intersects E (q) at 1 and q* « 1.768, 
and intersects L(q) at 1 and 1.761 (approx.). This line matters because the constant steady 
state is a strict local minimum for the energy when Bh X 2 < An 2 , by Theorem |10|(b). For 
example, suppose that 2 < q < 3 and h and X are such that E (q) < L(q) < Bh Q l X 2 < Aix 2 . 
Then the constant steady state h is a strict local minimum of the energy but is not a global 
minimum since £{h) > £{h ss ) by Theorem [ll](a) . A mountain pass idea then suggests the 
energy might have a saddle point at which its value is greater than £{h). Such a saddle 
ought to be an energy unstable positive periodic steady state, and should have period X 
and area hX. In fact we already know such a positive periodic steady state exists, by Rlffi 
Theorem 12] and Theorem |B| of this paper; we illustrate it in Figure [| Perturbing from the 
saddle in one direction leads to relaxation to the constant (see the top part of Figure §), while 
perturbing in the opposite direction yields apparent relaxation to a droplet (the bottom part 
of Figure |3|). The companion article fl9j contains more extensive numerical investigations. 



4. Proofs of Theorems |T| |5] 



4.1. Proof of Theorem [|. Recall that the positive periodic steady state h ss satisfies equa- 
tions (H and (El): 



(14) 



hZ + r(h ss )h' ss = 0, K s + H'(h ss ) = const. 



where H"(y) = r(y) is defined for y > 0. 

We compute the first four variations of the energy d). Let u G H\T X ) have mean value 
zero. For all e small enough that h ss + eu > (so that H(h ss + eu) makes sense), we have 



(15) 

(16) 

(17) 
(18) 



d 

de 

d 2 



S(h ss + eu] 



2 £{h ss + eu) 



de 



d 3 

—£(h ss + eu) 



6=0 



6=0 



6=0 



6=0 



[ [h'^ s + H'{h ss )]udx = 
Jo 

by (14), since u has mean value zero, 



x 



[(u 1 ) 2 - r(h ss )u 2 ] dx 



numerator of Rayleigh quotient ( ]TT| ) for Tx(h ss 
r'(h ss )u 3 dx, 



x 



r"(h ss )u dx. 



If the steady state h ss is linearly unstable with respect to zero-mean perturbations at 



period X, then the numerator of the Rayleigh quotient (111 ) is negative for some zero-mean 
trial function u G i/ 1 (Tx) \ {0}. We can assume u G C°°(Tx)- From (|I~6|), the second 
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variation of the energy in the direction u is negative, so that h ss is energy unstable in the 
direction u as desired. 



From Theorems 1 and 3 in [Tj| , h ss is linearly unstable with respect to zero-mean pertur- 
bations at period X if it is non-constant and: either the least period of h ss is X/j for some 
integer j > 2, or r = g/f is convex (r" > 0) and non-constant on the range of h ss . Hence h ss 
is energy unstable in those situations. 

Now assume h ss is non-constant and r is strongly convex (r" > 0). We consider the second 
variation of S in the direction u 



±Ks- 



&2 



S(h ss ±eh 







e=0 


Jo 1 



h 



ill 2 



r(h ss )h[ 



dx 



x 



r{h s 



dx by (14) 



x 



r(h ss ) [h'JC] dx 



x 



(19) 



x 



r'(h a 



h 



I 3 



dx = — 



r'(h ss )h' S sh" s dx 
x 

r"(h ss )h' ss dx. 



By assumption, r"(h ss (x)) > for all x. Since h' ss is not identically zero, the second variation 
in the direction u = ±h" s is negative, and so ±/i"s i s an energy unstable direction for h ss . 

It remains to prove h ss is energy unstable in the directions u = ±h' ss . Here the second 
variation is zero, since ([16]) becomes 



x 



[(K s ) 2 - r(h ss )(h'J 2 ] dx 



x 



[tC + r (^ss)^ s ] Ks dx = by parts and (gj. 



The third variation is zero for u = ±/igs ^y (|17|), because h ss is even about its minimum 
point while (±h' ss ) 3 is odd, by uniqueness for the nonlinear oscillator equation h" s + H'(h ss ) = 
const in (p!4|). The fourth variation is negative by (0), because r"(h ss (x)) > by assumption. 
Hence h ss is energy unstable in the directions u = ±h' ss , completing the proof. □ 

A number of other authors, working on closely related topics, have noted that translation 
invariance of the evolution implies the second variation of the energy in the direction h' ss is 
zero |], 1(| Those authors then proved energy instability in the direction u = h' ss + r] for 
some small function r]. Their arguments relied either on h ss having least period X/j for some 
j > 2 or else they did not impose a zero-mean requirement on the perturbation. Above, we 
have used instead the assumption r" > 0. 



4.2. Proof of Theorem g. 

When q < 1 or q > 2, Theorem follows from the last statement of Theorem [l] since 
r(y) = By q ~ l is strongly convex. 



20 



When q = 2 and u = —h" s , formula (|19| ) shows the second variation of the energy is zero, 
since r(y) = By and r" = 0. The third variation is negative because 



d 3 

j^£(h ss - eh'l s 



e=0 



B [ K) 3 dx by (0) 
Jo 

— 2£> / KsKsKs dx by parts 
Jo 

■2B- [ h ss (K s ) 2 K s dx since hZ = BhJ.L by (0) 
Jo 

~\& J ( h '^ dx <0 - 



Thus the steady state is energy unstable in the direction u = —h" s . 

Now suppose q > 1 and E'(a) > 0. To obtain the unstable direction u, start by rescaling 
k a to give 

P(a) 

(20) K a (x) := — tt rk a (P(a)x) . 

A[a) 

By construction, K a has period 1 and mean value 1. Then define 

d 

K a [x) := ——K a (x)] 
da 

K a is well-defined and smooth because P and A depend smoothly on a while k a (x) is jointly 
smooth in (a,x). Notice n a is even in x, has period 1, and has mean value zero: 

1 d f 1 d 

K a (x) dx = — — / K a {x) dx = — (1) = 0. 
da J oa 

(See [[18|, §5.4] for more properties of n a .) Let -u(x) = ±K a (x/X). 

We want to show u is an unstable direction. First we show that composing the rescalings 

of h ss to k a and then k a to K a yields 

(21) BX 2 h ss (xX) q - 1 = EK a {x) q ~ l . 

Indeed, from the definition Q2"0"D of K a , the righthand side of Q2"T| ) reduces to P 2 k a (Px) q ~ 1 , 
and then one can substitute for k a in terms of h ss using (|^). Next one obtains the lefthand 
side of ( pip by using the relation P = [Bj ' D) 1 l 2q D 1 / 2 X which relates the periods of k a and 
h ss (P and X, respectively). 

Now that we have (pl|) , use (0) to compute the second variation of the energy in the 
direction u(x) = ±K a (x/X) as 

f [n' a (x/X) 2 /X 2 -Bh ss {x) q - l n a {x/X) 2 ] dx = [ [^{xf - BX 2 h ss (xX) q - l n a (x) 2 ] dx/X 
Jo Jo 



l l [«) 2 - EK q - l K 2 a ] dx/X by © 
Jo 



< 
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by the proof of |L8|, Prop. 14], which uses both q > 1 and E'{a) > 0. Thus the steady state 
is energy unstable in the direction u = ±n a {x/X). □ 



4.3. Proof of Theorem |3|. Since E'(a) < by hypothesis, Proposition 15 of [18| implies 
> 0, where 



:= min 



fn \( v ') 2 ~ E(a)K^ l v 2 ] dx , f 1 

J ° — rl \\ a : v G H (Ti) \ {0}, / v(x)dx = 

J Q v z ax 



(In fact E\a) < would suffice to get this.) Notice = X 2 Ti(h ss ), by letting v(x) = 

u(xX) and using the identity (|2T|) and the definition (|TTD of T\. Hence Ti(h ss ) > 0. 

Consider u G i/ 1 (Tx) \ {0} with mean value zero. The first variation of £ in the direction 
u is zero by (|15l), and the second variation of £ in ([16]) is nonnegative because it equals the 
numerator of the Rayleigh quotient for Ti(h ss ). 

If the second variation is positive then h ss is energy stable in the direction h ss and we 
are done. If the second variation is zero then the Rayleigh quotient of u is zero, and so 



T~i(h ss ) = and u minimizes the Rayleigh quotient for T\ in (|iT|) . Hence u satisfies the 
Euler-Lagrange condition u" + r(h ss )u = const, and so it satisfies Cu = where C is the 
linearized operator defined in |1| eq. (4)] (take a = there). Theorem 10(a) in 0] and the 



hypothesis E'(a) < now imply u is a multiple of h' ss , and so u is odd. Therefore the third 
variation fllTf ) of £ in the direction u is zero, because u is odd and h ss is even. The fourth 
variation is positive by (0) because r(y) = By q ~ x and r"(y) < for 1 < q < 2. Thus h ss is 
energy stable in the direction u. □ 



4.4. Proof of Theorem |5|. As usual, the first variation of £ at h ss is zero by (|i~5|). We now 
prove non-negativity of the second variation, given by fli~6|). 

First note that h ss is symmetric about every point at which h' ss = 0, by uniqueness for 
the ODE h" B + H'(h ss ) = const (see (|T4|)); the uniqueness theorem applies here since the 
coefficient function H' is C l (even C 3 ) on the range of the positive bounded function h ss . 

Since h ss has a minimum at x = by hypothesis, we conclude h ss is even and h' ss > on 
(0, X/2) (otherwise h ss would have period less than X). 

Consider a minimizer of the Rayleigh quotient ([TT|) with respect to odd functions u\ it 
is a smooth odd function u satisfying u" + r(h ss )u + fu = for some constant f. Since 
u(0) = 0, one must have u'(0) ^ because otherwise u = by the uniqueness theorem for 
linear ODEs. Also, u(X/2) = by the oddness and periodicity of u, and so there is a point 
b G (0,X/2\ with u(b) = and ii / between and b. Assume u > between and b 
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(otherwise consider —u). Then 



r / uh' ss dx 
'o 



- / [u + r(h ss )u] h' ss dx since w" + r(h ss )u + tu = 
Jo 

-u'(b)h' ss (b) - / [ti£ + r(h ss )h'Judx by parts, since h' ss (0) = 
J o 



-u'C&KP) by (0) 



> 0, 



because u'(fo) < and h' SB (b) > 0. Since u and /ig S are positive on (0, b) it follows that f > 0. 
Hence the second variation of £ in ([Tj]) is nonnegative, as desired. 

The third variation at e = is zero by fll?!) , since u is odd and h ss is even. 



The fourth variation (18) is positive not just at e = but at all e G (0, 1): 
d 4 f x 

-— f(/i ss + eu) = — r"(h ss + eu)u 4 dx > 
cfe 4 Jo 

by the strong concavity of r, provided u ^ 0. Taylor's theorem completes the proof, since 
for some e G (0, 1), 

2! & 2 



S (h^ + u) = £ (hss) + ttt (^s + e«) 



e=0 



1 ^ <.„ 



> £(hss). 



□ 



5. Proofs of Theorems [6] — |TT| 
5.1. Proof of Theorem 19. We start by relating £(h ss ) to £(k a ). Take 



(22) 



H(y) 



»8+l 



2/ 



, 9^0,-1, 



9+1 

ylogy-y, q = 0, 
y-logy, q = -l, 



so that H"(y) = y q 1 , and recall from the definition (|J) of the energy that 



(23) £(/i s 



i.P(a) 

cfx and £(k a ) = 

Jo 



H{k a ) 



dx. 



Denote the period of h ss by P ss = X, and the area by A ss = J Q X h ss dx. Writing P = P 



a) 



and A = A(a), the rescaling ([7[) implies 

(24) 



P = l Kd> ss ' y ^ ' and 



e -B/2B B i/2p ss; g = 0, 
Notice that the rescaling (0) can be written as 



(%f 2q D^A ss , q^0, 

e -3D/2B B l/2 A ^ q = Q 
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From (|24 ) we obtain the invariance relation 

(26) BPj "At : = P 3 ^" 1 = E(a), 

and this implies 



(27) 



2 / p v 
p~ 



"A 



p 



Using (p5|) , (|26|) , (f27|) and the definitions (|23|) , we at last deduce a relation between £{h s 
and £{k a ): 



[2%) 



-1— ^ (^ ss ) ~ ^(^s) _ r~ -(9+1) £ (k a ) ~ £(k a ) 



p 



P(a) 



where the mean values are h ss := A ss / P ss and k a := A(pi)/P(a). (When checking (p8p, one 
can omit the linear terms in H(y) from the calculations, since h ss and /i ss have the same 
mean value, as do k a and k a .) 

In view of (p§|), then, Theorem || follows from: 

Proposition 12. Fix ai G (0,1). 

If q > 2 or q < 1 then £(k a ) > £(k a ) for all a G (0, 1). 
If q = 1 then £(k a ) = £{k a ) for all a G (0, 1). 

Ifl<q<2 and E'(a) > Va G (a x , 1) then £{k a ) > £(k a ) Va G [«i, 1). 
Ifl<q<2and E'(a) < Va G (c*i, 1) t/*en £(jfe a ) < £{K) Va G [«i, 1). 

Note that if g > 2 or g < 1 then E'(a) > for all a by ||T8| , Theorem 11]. 



Proof of Proposition [71. The proof depends on a number of elementary differential equations 



and inequalities that we derived in §§5.1-5.2 of flafl , and the reader may wish to skim those 
sections before proceeding. 

If q = 1 then £(k a ) = £(k a ), as one sees directly from the formula in (|23|), using that 
k a (x) = 1 + (a — 1) cosx, P(a) = 2n and k a = 1. So we assume g ^ 1 from now on, and 
a G (0,1). 

By definition, 



£(k a ) £{ka) 
P 

(29) 

First assume q ^ — 1; then 

d — £(/Cq,) 



1 

P 



^« - H(k a ) 



dx + H{A/P) 
±£ k' a 2 dx-H(a) + H(A/P) by [18, eq. (21)]. 



(30) 



da 



P'(a 

p2 



£ K 2 dx + (l + (q + l)H(A/P) 1 j (± 
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where we have used |H| eq. (35)] and the identity H'(y) = 1 
q —1). Differentiating the function 



l)H(y)y 1 (valid for 



(31) 



^(«) : = ( J 



"(9+1) 



£{ka) £{k 
P 



(which is inspired by (|2"8"D) with respect to a, we find from fl29p and (^) that 



.4 



P 

Substituting 



"(9+2) 



1 



"(9 + 1) ('^ 



AP' 

~p2 



£ + (q + l)H{a) 



A , -t \ t t / \ 9 + 31 
_ + (g+l)P( a ) = l__ 



from []T8|, eqs. (31-32)] yields 



(32) J* (a) 



1 (A 



2 VP 



"(9+2) 



1 



(9 



££ 2 da; 



^ ) + 2.4 P-" 2 P' 



for q ^ — 1. 



For g = -1 we obtain exactly the same formula (|32|) for J 7 ' (a), as follows: 
d £(k a ) — £(k a ) 



da 




p 


P' 

~P2 


i: 








1 1 






p 



pi fP /A V 

— y A£ 2 da: + (1 - P/A) f-j from (H) and 0, eq. 



(35)] 



k' 2 dx 



-2 



+ 2AP~ 2 P' 



using that (1 — P/A) = J Q (k' a ) 2 dx/A when q = — 1, by JIB], eq. (28)]. The last equation is 
(0) for g = -1. 

Equation ( |3"2"D simplifies to 

~1_ , / A x ~ 29 ^ 



(33) 



-P~ 



P 



fe/-v Q/00 



E'(a) 



using that P = P^A/P) 9 " 1 by fl2§). Hence 

(34) E' > P' < and P' < J"' > 0. 
Also 

(35) J-'(a) — > as a — > 1, 

by the formula (|29D together with the facts that P, v4 — > 27r and P/A — > 1 as a — > 1 (see 



fT8| , Lemma 6]) and that k' a — > uniformly as a — > 1, by [18, eq. (21)]. Proposition [12] now 
follows from fl34|) , (j35[) and |18|, Theorem 11] (which shows P' > when q > 2 or q < 1). 
For example, if P' > on (ai, 1) then T' < on (ai, 1); since P(l) = we deduce T > 
on [ai, 1), and so £(/c Q ) — £{k a ) > for a e [«i, 1). □ 
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5.2. Proof of Theorem [7|. The proof relies on formulas for A' (a) and E'(a) that were 
derived in Lemmas 16 and 18 of for all q ^ — 1, 



(36) 



A' = -(g+l)ff(a)P' 



-H'(a)P, 



(37) 



9-2 



P' [(g - 3)A + (g - l)(g + l)H{a)P] + l -{q- l) 2 H'(a)P 2 



Note that E' > if q < 1 or g > 2, by [TTSL Theorem 11]. Our assumptions therefore imply 



-1 < g < 3 and E'(a) > for all a £ (0, 1). Also £(0) > because g > -1 0, §3.1.2]. If 
we define 



X = [E(0)/(BAI~ 



g-lAl V(3-9) 



then 



BX 3 " 9 ^- 1 = P(0) < S(a) = BX 3 " 9 ^- 1 

by (pB[); here the value a £ (0, 1) is determined by translating and rescaling h ss to a particular 
k a , as in § |2~3| . Hence < X < X (using that q < 3). 

By rescaling the zero contact angle function ko on the interval [0, -P(O)] (as in §[T3 



sec 



|T7| , Claim 5.1.2] for details) we obtain a zero contact angle droplet steady state h ss of 
h t = —{h n h xxx ) x — B(h m h x ) x , with length X and area A ss as desired. 

It remains to prove that the energy of this droplet steady state is lower than the energy 
of the positive periodic steady state h ss . That is, we want to prove 



x r 



\{K s f 



BG(h s 



dx < 



x 



where 



»4+l 



BG(h s 



dx 



G(y) 



9(9+1) ) 

y\ogy-y, q = 0, 
-logy, g = -l; 



note that we can use G instead of H in the energy because they differ only by a linear 
function (cf. the definition (|22"D of H) and h ss and h ss have the same area, A ss . 
Since h ss is supported on (0,X) and because G(0) = 0, the desired inequality is 



(38) 



x r 



(K s y - BG(h s 



dx < 



x r 



/ \2 



BG(h s 



dx. 



Next rescale h ss to k and h ss to k a : from and ( |27| ) (with P ss replaced by X or X as 
appropriate) we deduce ( |38| ) is equivalent to 



-4, 



A(0)J 



P(0) 
X 



P(0) 



dx < 



A R 



A(a) 



P(a) 
X 



P(a) 



/ \2 



G(k c 



dx, 
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except that when q = we have to subtract 



log 



a ss p(oy 

A(0) X . 



kn and 



log 



^4 SS Pja) 
A(a) X 



k n 



from the integrands on the left and right sides, respectively. By substituting the relations 
X =[E(0)/(BA g - 1 )] 1/{3 ~ q) and X = [E{a)/{BA£ X ^ 1/(3_,) 



into the last inequality and using the definition E(a) = P(a) 3 q A(a) q , we see the desired 
inequality reduces to Q(0) < G(ot) where (for q ^ 3) 

»P(a) 



(39) G(a) = A(a 



,(9+3)7(9-3) 



(KY - G(k a ) 



dx 



log A(a), when q = 



Thus to show the energy of h ss is lower than that of h SB , it suffices to show Q'{ot) > for all 
a G (0, 1), assuming — 1 < q < 3 and E'[a) > for all a G (0, 1). 



To show Q' > 0, we substitute [18|, eq. (32)] and [[L8|, eq. (30)] into the definition fl39|) of 
Q, obtaining that 



when q ^ 3, 0, —1, —3. 



q(q + 3)" ' q + 3' 

After differentiating this formula with respect to a and then substituting for A' (a) from 
(P^|), we simplify with the help of ( |3"7| ) to obtain 



(40) 



& = J-H ia)A */i*-*) (± 
q-3 \P 



-q+2 



E' when q=£ 3,0, -1,-3. 



When q = we obtain (/ = | — ^HPA 1 + | log A, by putting the q = versions of [18 
eq. (32)] and @, eq. (30)] into ©. Hence ^' = -\HP- 2 E', by differentiating and using 
(H) and (|3). That is, (0) holds when q = 0, also. 



We conclude from ( [40]) that Q 1 > as desired (provided — 1 < g < 3 and E'(a) > for all 
a G (0, 1)), since for q > — 1 we have H(a) < by the definition (p2|) . □ 



5.3. Proof of Theorem The proof involves rescaling arguments rather than the energy. 
Write P ss = X. Assume h ss is non-constant. 

Suppose that in fact h ss is positive (and so smooth by hypothesis). Then the least period 
of h ss equals P ss /j for some positive integer j, with the area per period equaling A ss /j. If 
j = 1 then h ss must be a translate of h ss , by modifying slightly the uniqueness remarks in 
flB] , §6.2] (using the assumption that E' < to get strict monotonicity of E). 

Thus we can assume j > 2. By rescaling h ss and h ss to fc a and fca for some a, a G (0, 1), 
as in (0), we get from fl2"6"D that 



E(a) = BPt q A q s ; 1 and E(a) = fi(P ss /j) 3 - g (A ss /j) 



9 -l 
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Hence 



2 E(a) £ 
A<j 2 = -^-{ < since £' < 

~ J E(a) E(l) 



12, , / 1 1\ 3 9 / 3 l x<? 1 



l + — ,- £? — ,- =:J(g)say, 



4vr 2 g ^ 7 \ 2 q 2 J W 2 
by the formula for £(0) in 0, eq. (3.13)] and since £7(1) = P(1) 3 -*A(1) ?_1 = 4vr 2 by p 



Lemma 6]. We will obtain a contradiction by showing J(q) < 4, when 1 < q < 1.75; this 
will show h ss is not positive. 
For 1 < q < 1.5 we have 



12, . / 1 1\ 3 1 / 3 l xl " 1 



^(9) < — + 1-5) -B , - B[ ,- ^3.17<4, 

w; -4vr 2 r ; V2 - 1.5 2 y \2-1.5 2/ 

where we have used that the Beta function B(a, b) = t a_1 (l — t) b ~ l dt is decreasing in its 
arguments, and is bigger than 1 when those arguments are less than 1. For 1.5 < q < 1.75 
we similarly have 

12, _ , J 1 xx 3 - 1 - 5 / 3 



J(q) < — (1 + 1.75)5 , - B\ , - «1.73<4, 

W; -4vr 2 1.5 V ; V 2 ' 1 - 75 2 J V 2 " 1 - 75 2 J 

completing the contradiction. (Actually the best upper bound for J(q) seems numerically 
to be about 1.04, when 1 < q < 1.75; see [|TJ|, Figure 6].) 

The above contradiction implies h ss is not positive everywhere. Consider therefore one 
component of the set {x : h ss (x) > 0}, say an interval with length P < P ss . Write A < A ss 
for the area under h ss on this interval. Note the contact angles of h ss must be the same at 
the two endpoints of the interval, as a consequence of the nonlinear oscillator equation (§) 
(see for example JL7|, §2.2]). Suppose these contact angles are zero, so that h ss rescales to k$ 
on the interval, using (]?]). Then 

£(0) = BP z - q A q - 1 < BP*- q A q - x = E{a) 

by applying (p6|) twice, but this contradicts our assumption that E' < 0. Thus the contact 
angles of h ss must all be nonzero, as desired. □ 

5.4. Proof of Theorem |9|. Translate h BB x and h ss2 so that they attain their minimum values 



at x = 0, and then rescale as in § |2.3j to obtain k ai and k a2 respectively. Since h ss i and h ss2 



have the same period and area, for which we write P ss = X and A ss respectively, it follows 



from flUD that £(«i) = E(a 2 ). Notice h ss \ 7^ h ss2 ^> «i / a 2 , in view of the expression 
(f2~5|) for h ss in terms of k a and P ss , A ss , P(a), A(a). Since E(a±) = E(a 2 ) while £ is strictly 
decreasing on (0, a. C rit) and strictly increasing on (a cr i t , 1), we conclude a cr u must lie between 
«i and a 2 . 
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We show a\ < a 2 - The hypothesis min/i ss i < min/i ss2 gives 

D 1 \ 1/q ^ (D 2 ^ lq 



by the rescaling (|7|). Next apply the first equation in ( p4| ) to solve for Di in terms of P(ai), P s 
and 23, and similarly solve for D 2 in terms of -P(«2), -Rs and B. Substituting into the above 
inequality gives aiP(ai) 2 ^ q ~ 1 ^ < a 2 P(a 2 ) 2 ^ q ^ 1 \ The strict increase of a \— > aP(a) 2 ^ q ^ 1 ^ 
implies a,\ < a 2 . 

Since a\ < a crit < a 2 , our hypothesis on E' implies E'(a\) < and E'(a 2 ) > 0. Theorem || 
then implies that h SB x is energy stable, and Theorem |2| implies h ss2 is energy unstable. 

Next we show S(h ssl ) < £(h ss2 ), or £(h ssl ) - £(h ssl ) < S(h ss2 ) - S(h ss2 ). In view of the 
rescaling relation (p8|) for the energy, it suffices to prove F(ai) < jF(a 2 ), where T was 
defined in (j3l"l). 

To prove T(a{) < F(a 2 ), we write (^) as J 7 ' (a) = 5(a) (l/E)' (a), where 

S3/1-2 u 2 . 



5(a) = —P A I k'J dx. 







The point of this transformation is that 5(a) is strictly decreasing: P' < and A 1 > 0, by 
|TT| , Props. 7.3 and 7.4], while a i— > J Q P fc^ 2 is decreasing by [18|, eq. (35)]. Also (l/E)' > 



on (0, otcrit) and (l/E)' < on (a cr it, 1), by assumption. Hence 

F'(a) > 5(a crit )(l/E)'(a) for a G (0, l),a 7^ a crit . 

Integrating this inequality from a± to a 2 and using that E(a\) = E(a 2 ) gives J : (a 2 ) > F(ai) 
as desired. 

Finally, © shows ^'(a) < on (a crit , 1), and so .F(a 2 ) > T(l) = by fl35|). Thus (|8| 



yields £ ( h ss2 ) > £ ( /i ss2 ) . □ 



5.5. Proof of Theorem [T0| . From the definition ( IT ) we see 



Ti(n) = mm — — r(/i) = I — ) — r(n), 



J Q X it 2 dx 

with the minimum being attained precisely for linear combinations of sin(27rx/X) and 
cos(2kx/ X). The first two paragraphs of the theorem follow directly. 

(a) Now suppose r(h)X 2 > An 2 . The variational formulas (|TJ|) and ( |16D in the proof 
of Theorem [l] show the constant steady state h is energy unstable in the sine and cosine 
directions, since these are Ti(/i)-eigenfunctions and Ti(h) < 0. Suppose next r(h)X 2 = Att 2 
and r"(h) > 0. Then the first two variations of the energy in the ±sin directions are zero, 
by (|TJ) and (|i~6|) . The third variation equals r'(h) times the integral of =Fsin 3 , by (|17|) ; thus 



the third variation is also zero. The fourth variation is —r"(h) J Q X sin 4 (27rx/X) dx, which is 
negative because we assumed r"(h) > 0. Thus the constant steady state is energy unstable 
in the ± sin directions. Argue similarly for the ± cos directions. 
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(b) If r(h)X 2 < 47r 2 , or if r(h)X 2 = An 2 and r"(h) < 0, then n(ft) > and so we get 
energy stability by modifying the argument of part (a) as follows. The first variation of the 
energy in a direction u is always zero. If the second variation is positive then we are done. 
Otherwise the second variation must be zero, so that T\{K) = and r"(h) < 0. Then u must 
be a linear combination of sines and cosines, and so the third variation is also zero. Then 
the fourth variation is positive. 

Furthermore, for r(h)X 2 < 4tt 2 we will prove h is a strict local minimum of the energy, 
and is nonlinearly stable in H 1 . In doing this, we will use below a certain sufficiently small 
number 5 £ (0, 1). Then considering u £ H l (Tx) with mean value zero and ||w||iyi 



we find for all e E [0,5] that 

j£(h + eu) = / [(u') 2 -r(h + eu)u 2 ] dx 
Jo 



d 2 
de 2 ' 



> 5\\u 



5 [ [{u'f + u 2 ] dx + {l-5) [ 
Jo Jo 

S) 



\H 1 (T X ) 



+ 1 



47T 



f\2 



r(h + eu) + 5 



X 2 



u 



1-5 

2 



-U 



dx 



dx > S\\u\\ H 



1 (Tx)' 



In the second-to-last step above, we used r(h)X 2 < An 2 and ||w||oo < Cll^ll-f/ 1 and we chose 
5 sufficiently small (independent of u and e). 

On the other hand, \eu\ < 5\\u\\l°° < h/2 provided 5 is chosen small enough (independent 
of it), and hence 



d 2 
de 



2 £(h + eu) 



x 



[(u'f - r(h + eu)u 2 ] dx < C(h) / [{u'f + u 2 ] dx 

Jo 

for some constant C(h) > 1 
= C(h)\\u\\ 2 Hl{Tx) . 



We deduce from the preceding estimates and the vanishing of the first variation in ( |i~5D 
that if u £ H l {Tx) has mean value zero and ||m||#i(t x ) < 5, then h + u > and the energy 
varies quadratically away from h: 



(41) 



<e(h+u)-e(h)<-c(h)\\u\\ 2 HHTx) . 



The lefthand estimate implies h is a strict local minimum of the energy, with respect to 
X-periodic zero-mean perturbations. 

We now prove h is nonlinearly stable, in the sense that if h(x, t) is a smooth positive 
solution of ([!]) for x £ Tx and t £ [0, T], for some T > 0, and if h(-, 0) has mean value h and 
||/i(-,0) - h\\ H i(T x ) < Js 3 /4C(h)i then \\h(-, t) - h\\ HHTx) < 5/2 for all t £ [0,T]. Indeed, 
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the quadratic bounds (§TJ) and the dissipation of the energy together imply 

h\\h(;t)-h\\l 1{Tx) < £(h(;t))-£(h) 



< £(h(;0))-£(h)<lc(h)\\h(;0) 



ni2 153 

n \\m(T x ) < 



so that ||/i(-, t) — /i||fl-i(T x ) < 5/2 for all t G [0, T}. This stability result holds for all sufficiently 
small 5. □ 

This proof of nonlinear stability for a linearly stable constant steady state does not carry 
over to linearly stable positive periodic steady states h ss , because there Ti(h ss ) = (due to 
translational null directions). This zero eigenvalue is absent for the constant steady state, 
since translation of h gives h again — a trivial perturbation. Imposing Neumann boundary 
conditions eliminates the translational perturbations and their associated zero eigenvalue. 
Hence, the nonlinear stability proof would hold for a positive steady state under the Neumann 
boundary conditions h x = h xxx = at x = 0,X/2 (see §|2.6|), provided the steady state is 
strictly linearly stable, i.e. the first eigenvalue is positive. 

5.6. Proof of Theorem [TT| . Suppose first that — 1 < q < 3. By (|26| ) with a = 0, a steady 
state h ss supported on a single interval of length X < X and with area A ss = hX and zero 
contact angles exists if and only if for some length X < X we have 

(42) BX z - q Al; 1 = E{0) = P(0) 3 - ? A(0) <? - 1 . 

That is, if and only if BX z - q (hX) q - 1 > E(0), which is (0). 

Suppose such a steady state h ss exists, supported say on the interval (0,X). By above, 

X = B l/{q - z) P(Q)[A(<d)/hX}^l^ q \ 

We want to show that £(h ss ) < £{h) if and only if (|13D holds. For this, compute using G 



like in § |5.2j to find 

£(^ss) 



£{h) 



x 



x 



hi 



BG(h ss ) + BG(h) 



dx 



(43) 



P(0) 



BG(h ss ) 



dx + BG{h)X 



G(k 



dx + BG{h)X 



A(0) 2 X 3 

when q > — l,q ^ 0, by rescaling from h ss to k and using (|25|) and (|27|) (with P ss replaced 
by X, and A ss replaced by A ss ). Putting a = into |L8|, eq. (32)] and eq. (30)] and using 
H(0) = enables us to evaluate f (k' ) 2 dx and J G(k ) dx, and hence we deduce 

i s 2 s P(0) 3 q-3 



(44) £(h ss )-£(h) 



A(0) 2 X 3 q{q + 3) 



A(0) + BG(h)X when q> -l,q^ 0. 
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After substituting the definitions of X and A ss = hX from above, we find 



£(h ss ) - £(h) = -B^^- q \hX)^l^ 

q 



U(o) 2 «/fo- a >+— —(Bh q ~ i x 2 y 



_g + 3 q + 1 

when — 1 < q < 3, q ^ 0. Plainly now S(h ss ) < S(h) if and only if ( |l"3"l) holds. 
When q = we find has an extra term -BA ss \og A ss P(0)/A(Q)X 



so that after 



again using |18], eq. (32)] and [18|, eq. (30)] and substituting for X and A SSl we obtain 



- 1 — A(0) 2 /e 
S(h ss ) -£(h) = -BhX log __{ 1 when o = 0. 

3 Bh X 2 



Remembering that A(0) = 2e 3//2 A/vr/3 when g = 0, from |[7|, §3.1.2], we conclude £(h ss ) < 
£(h) if and only if Bh 1 X 2 > 4e 2 ir/3. Incidentally, one can check L(q) is continuous at 
q = 0. 

For g > 3, simply modify the above proof from the case —1 < q < 3. Notice when q = 3 
that (|I2"D becomes B(hX) 2 = E(0), which yields no formula for X. And when g > 3 we get 
S(h ss ) > £(h) because the first term in ([44]) is nonnegative and the second is positive. □ 

6. Conclusions and Future Directions 

If you perturb a positive periodic (or constant) steady state h ss of the evolution equa- 
tion ([!]), without changing its area, then towards which steady states might the solution 
subsequently evolve? 

That is the broad question addressed by this paper. To answer it, we focused on three 
specific questions: 

• Existence: Do there exist other steady states having the same area and same period 
as h ss , or having period a fraction of the period of h ss 7 If so, then these other steady 
states are plausible contenders for the long time limit. The constant steady state h ss 
obviously satisfies these requirements, but there might be another positive periodic 
steady state (different from h ss ) that does also, or perhaps a array of droplet steady 
states having zero or nonzero contact angles. 

• Stability: Are h ss or these other steady states linearly stable? energy stable? If a 
steady state is the long time limit of some generic solution, then surely it must be stable 
under perturbations. 

• Relative energy levels: Do any of these other steady states have lower energy 
than h ss 7 Only steady states with lower energy are accessible, when starting from a 
small perturbation of h ss . 

The existence question was substantially answered for power law coefficients by |18, The- 



orem 12], Theorems [7] and [| and [Tj| Figures 3-5], also using Theorem [Tl| when h ss is the 
constant steady state. But the existence question remains open for droplet steady states 



32 



with nonzero contact angles, if we wish to specify the area and the length of the support. 



Some information on nonzero angle droplet steady states is in our earlier paper [[171 , §5.2]. 

The stability question was resolved for positive periodic steady states in the power law 
case by fl8| , Theorems 1,3,7,9] and Theorems [l]-[5] here. In particular, Theorem [I] related 
linear instability to energy instability. Theorem [1^ handled the case of constant steady 
states. Our numerical simulations in the companion article |19| suggest that linearly unstable 
steady states are indeed nonlinearly unstable, with the linear behavior dominating near the 
steady state, but this observation is so far unsupported by a 'linearization theorem' for the 
power law evolution h t = —{h n h xxx ) x — B(h m h x ) x . (A linearization theorem is known in the 
Cahn-Hilliard case / = 1, by using semilinear operator theory; see for example §6].) 

The energy level question has been largely settled in the power law case by Theorems [| 
|7J and |S] when h ss is positive and periodic, and by Theorems |5] and |TT] when h ss is constant. 
When h ss has compact support with zero contact angle, use Theorems ^ and [II]. For example, 
Theorems |9| suggests that when m = n + 0.77, there can exist two positive periodic steady 
states with the same period and area, and that the unstable steady state has higher energy 
than the stable one. Our numerical simulations have found heteroclinic connections from 
the high energy steady state to the low energy one. 

Future directions. The stability question for droplet steady states (with zero and nonzero 
contact angles) is open. So is the problem of computing relative energy levels of non-zero 
angle droplet steady states with respect to zero-angle droplets and constant and periodic 
steady states. 

Another open problem is to answer the existence, stability and energy level questions for 
general coefficient functions / and g. We have treated power law coefficients, and Grinfeld 
and Novick-Cohen [14| cover the Cahn-Hilliard equation, with / = 1 and g(y) = 1 — 3y 2 . 
But for general coefficients, about all we know is that every positive periodic steady state is 
linearly and energy unstable when g/f is a convex function, by Theorem [I]. 

Finally, many of our existence, stability and relative energy level theorems for the power 
law evolution would be improved if we knew E'(a) < for all a when 1 < q < 1.75. We have 
not been able to prove this conjecture, though numerically it is clear from [IR Figure 3]. 
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